Quantum master equation for the microcanonical ensemble 
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By using projection superoperators, we present a new derivation of the quantum master equation 
first obtained by the Authors in Phys. Rev. E 68, 066ff2 (2003). We show that this equation de- 
scribes the dynamics of a subsystem weakly interacting with an environment of finite heat capacity 
and initially described by a microcanonical distribution. After applying the rotating wave approx- 
imation to the equation, we show that the subsystem dynamics preserves the energy of the total 
system (subsystem plus environment) and tends towards an equilibrium state which corresponds to 
equipartition inside the energy shell of the total system. For infinite heat capacity environments, 
this equation reduces to the Redfield master equation for a subsystem interacting with a thermostat. 
These results should be of particular interest to describe relaxation and decoherence in nanosystems 
where the environment can have a finite number of degrees of freedom and the equivalence between 
the microcanonical and the canonical ensembles is thus not always guaranteed. 

PACS numbers: 05.30.-d; 03.65.Yz; 76.20.+q. 
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INTRODUCTION 



Irreversible processes do not occur in isolated quantum 
systems with a finite number of quantum levels. In order 
to understand relaxation toward equilibrium in these sys- 
tems, one needs to take into account the effect of their in- 
teraction with their environment. The usual way to pro- 
ceed is to consider a total Hamiltonian system composed 
of a subsystem coupled to an environment. The subsys- 
tem dynamics is then described by the reduced density 
matrix of the subsystem which is obtained by tracing out 
the degrees of freedom of the environment from the den- 
sity matrix of the total system. In order to drive the 
subsystem into an effective irreversible process over rea- 
sonably long-time scales (of the order of the Heisenberg 
time scale of the environment) on need to assume a quasi- 
continuous environment spectrum. This is typically valid 
when the energy spacing between the unperturbed quan- 
tum levels of the total system which are coupled together 
by the subsystem-environment interaction is sufficiently 
small to make the interaction between these levels effec- 
tive enough to "mix" them [2(. Once this condition is 
satisfied, the generic way to obtain a closed master equa- 
tion for the reduced density matrix of the subsystem is 
to use perturbation theory and the Born-Markov approx- 
imation which implicitly rely on the assumption that the 
environment has an infinite heat capacity and cannot be 
affected by the system dynamics. The environment thus 
plays the role of a thermostat for the subsystem and the 
descriptions in the canonical and microcanonical ensem- 
bles are equivalent Q. The resulting quantum master 
equation is called the Redfield equation 0, H, H, 0; Hi- 
Since the Redfield equation can break the positivity of 
the subsystem density matrix it is sometimes simplified 
further by time averaging the equation in the interaction 
representation (rotating wave approximation) in order to 
get a master equation of Lindblad form which preserves 

positivity @, m no, m m . 



In the present paper, we consider situations where the 
environment has a quasi-continuous spectrum but a fi- 
nite heat capacity. This means that the energy quanta of 
the subsystem may significantly affect the microcanonical 
temperature of the environment so that the equivalence 
between the microcanonical and canonical statistical en- 
sembles is compromised and only the microcanonical en- 
semble can be used a priori. Such situations should be 
generic since the density of state of a system growth expo- 
nentially with the number of degrees of freedom whereas 
the heat capacity only growth linearly. This means that 
there exists a range in the number of degrees of freedom 
where the quasi-continuous assumption can be satisfied 
without necessarily implying an infinite heat capacity. 
This domain where kinetic processes can occur in the 
subsystem but for which the usual master equations fail 
should be of particular importance in nanosystems where 
the number of degrees of freedom constituting the envi- 
ronment is not always large enough to be supposed infi- 
nite. In this sense, the present work can be viewed as a 
contribution to the recent att emp ts t o app ly statistical 
physics to small systems [3 El El EJ EE El • Our 
results should also be of relevance for systems which dis- 
play negative heat capacities in the thermodynamic limit 
and to which the microcanonical ensemble description 
applies. A well known example is provided by coupled 
spin systems because their spectra remain bounded and 
forms bands of finite extension [H HH HI H3- Other- 
known examples of such systems are reported in Refs. 

[II [IE [II [la. 



The plan of the paper is the following. In Section[TTl we 
give a new and more enlightening derivation of the equa- 
tion of Ref. [l[ by using projection superoperators. In 
Section IITT1 we apply the rotating wave approximation to 
our equation and prove that it conserves the energy of the 
total system. We also demonstrate that the subsystem 
relaxes to an equilibrium distribution which corresponds 
to equipartition inside the energy shell of the total system 
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as expected in the microcanonical ensemble. In Section 
IIV[ we show that, for large heat capacity, our equation 
reduces to the Redfield equation and we discuss in detail 
the conditions under which the master equations for in- 
finite thermostat fail to describe the correct subsystem 
dynamics and need to be replaced by our master equa- 
tion. Conclusions are drawn in Section fVl 



II. NEW DERIVATION 

In this section, with the help of projection superopera- 
tors, we give a different derivation of the quantum master 
equation first derived in Ref. 

We consider the dynamics of a system with a Hamil- 
tonian of the form 



H = H + XV , 



(1) 



where A is a small dimcnsionless parameter. The density 
matrix p(t) of this total system obeys the von Neumann 
equation 

P {t) = Cp(t)=~[H,p(t)] (2) 

= (C a + XC')p(t) = - l -[H ai p(t)}-X^[V, P (t)} . 

The solution of the von Neumann equation reads 

p(t) = W(i)p(0) = e ct p(0) (3) 

In the interaction representation where 

Pl (t) = e- Cot p(t) = e* Hot p(t)e-t Hot , (4) 
C'j(t) = e- Cot C'{t)e Cot , (5) 

the von Neumann equation takes the simple form 

Pl (t) = XC' t (t) Pl (t) = -\\[Vi{t), Pl {t)] . (6) 

By integrating Eq. ^ and truncating it to order A 2 , we 
get the perturbative expansion 

Pl (t) = W(t)p(0)=e- c " t e ct p(0) (7) 
= [W (t) + XWi(t) + \ 2 W 2 {t) + 0(X 3 )] p(0) , 



where 



Wo{t) = i; 



(8) 



dT£'j(T) 



W 2 (t) = ( dT [ dr L'j(T)L'j{T - t) . 
Jo Ja 

The inverse of W(t) reads 

W^it) = W (*)-AWi(i) (9) 
+A 2 [Wl(t)-W 2 {t)} +0(X 3 ) . 



Indeed, one can check that W^W -1 ^) 
For later purpose, we also notice that 

VV(t)-4W _1 (t) = XWi(t)A 



1 + G(A 3 ). 



(10) 



+A 2 



W 2 (t)A - VVi {t)AWi (t) + O ( A 3 ) 



Now, we consider a subsystem S interacting with its 
environment B. The Hamiltonian of this system is given 
by Eq. Jl} where 



Hn — Hq 



(11) 



S K (resp. B K ) is a coupling operator of system S (resp. 
B). We will use the index s (respectively b) to label the 
eigenstates of the Hamiltonian of system S (resp. B). 

We will use Liouville space where operators are 
map ped into vectors and superoperators into matrices 
[28j |. We recall some basic definitions 

scalar product : ((A\B)) = tvA^B , (12) 

identity : X = ^ \nn'))((nn'\ , (13) 

\nn')) «-> \n)(n'\ , ((nn'\ *-> \n')(n\ . (14) 

Useful consequences of these definitions are 

((nn'\nn')} = S nn 8 n >n' (15) 
((nn'\A)) = (n\A\n') . (16) 

An operator A in the total space which can be written 
as a product of a system and reservoir operator A — 
As® Ab will read in Liouville space \Aj) = \As)) <8> |Ab))- 
We now define the following projection superoperators 
which act on the Liouville space of system B 



V 



E 



Q = ^2\bb'Ubb'\(l-6 bb/ ) 



(17) 
(18) 



b.b' 



and which satisfy the usual properties of projection su- 
peroperators V + Q = 1b,V 2 = V,Q 2 = Q and 
VQ = QP = 0. V, when acting on the density matrix of 
the total system, eliminates the environment coherences 
but keeps the environment populations and leaves unaf- 
fected the subsystem degrees of freedom. Similar pro- 
jection superoperators have been recently used in Refs. 
[3d l3ll. I32I [33l] . For comparison, the projection superop- 
erator used to derive the Redfield equation reads 



VRed — 



b 



P 



(19) 



where p eq is the equilibrium density matrix of the envi- 
ronment. The two projection superoperators act differ- 
ently on the density matrix p(t). In the Liouville space 
of the total system, we get respectively 



v\ P (t))) 

VRedW))) 



£«M>|p(*)» c 

b 

\Ps(t))) ® \ P eq )) 



(20) 
(21) 
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((bb\p(t))) is a density matrix in the system space which 
depends on the environment state b. V therefore cor- 
relates the system state with the environment state. On 
the contrary, pR e d assumes that the system reduced den- 
sity matrix ps{t) = J2 b ((bb\p{t))) = trsp(i) is indepen- 
dent from the environment state which always remain at 
equilibrium. 

We now let V and Q act on the density matrix of the 
total system in the interaction form ([7]) and find 



V\ P i{t))) 

Q\pi(t))) 



VW(t)(P- 
QW(t)(V 



Q)\pM)) 

C)|pi(0)» 



(22) 
(23) 



From now on, we will consider initial conditions such that 
Q\p(0))) = 0. This means that the environment part 
of the initial condition is diagonal in the environment 
eigenbasis and is thus invariant under the evolution when 
A = 0. Taking the time derivative of Eq. (|2"2"|) and Eq. 
HMD and using |p/(0))) = W" 1 (t) | pi (i))), we get 

P\f>i(t))) = VW{t)VW-\t)V\ Pl {t))) (24) 
+VW{t)VW- 1 (t)Q\ Pl (t))) 

Q\pi(t))) = QM^VW-'^Vlpjit))) (25) 
+QW{t)VW- 1 {t)Q\ Pl {t))) . 

These equations are still exact. If we restrict ourselves 
to second-order perturbation theory, we can obtain the 
important result that the V projected density matrix evo- 
lution is decoupled from the Q projected part. Indeed, 
with the help of Eq. (fTU|) , we have 

VW{t)VW- x {t)Q = XPWi(t)VQ (26) 
+\ 2 VW 2 (t)VQ - X 2 PW 1 (t)PW 1 (t)Q + C(A 3 ) . 

The two first terms of the right-hand side are zero be- 
cause PQ = and the third one also because 



vm(t)r = ^{bbHbbic'Mb'b'Ub'b'i 



(27) 



b.iy 



= 4E 5K wE iwi wu iwi . 

k b.b' 

where (6| [B K {t), \b')(b'\] \b) = 0. 

After having showed that the relevant projected den- 
sity matrix evolves in an autonomous way, we will now 
evaluate the generator of its evolution using second-order 
perturbation theory. Again using Eq. (|10[) , we find that 

PW(t)PW^(t)? = XVWi(t)V (28) 
+\ 2 VW 2 (t)V - A 2 -pWi(i)7>Wi(i)P + C(A 3 ) . 

The only term of right-hand side which is not zero is the 
second one [see Eq. ([2"T])] whereupon we get 



V\pM)) = A 2 P / drC'MC'^t-TyP^t))) 
Jo 

+0{\ 3 ) . 



(29) 



Now leaving the interaction representation and using the 
fact that Ve~ Cot = e~ Cst P, we obtain 



T\ P (t))) = C s V\p{t))) 



(30) 



+X 2 e Lst P / drL'^C'jit - r)e- Lst V\p{t))) + 0(X S ) . 
Jo 

Now, wc define the quantity P(E b ,t) that will become 
the fundamental quantity of our theory 



P(E b ,t) 
n{E b ) 



((bb\r\p(t))) 



(31) 



where n(E b ) = txB${E b — Hb) is the density of state of 
the environment. Equation (|30|) can thus be rewritten as 



(32) 



P(Et,t) _ Cs P{E b ,t) , n2 



n{E b ) 



n{E b ) 



h , Jo 



4(tKK*-^)l^> 



n(E b .) ' 



where we stop explicitly writing +0(A 3 ) from now on. 
Equation (f3"2")) can be evaluated further and we get 



P{E b ,t) = P(E b ,t) A 
n(E b ) S n(E b ) //-'— — , 



(33) 



■****<&! [s K (t)B K (t), [s K '(t - T )B K , (t - t), 

t |6 / ){&'|ll|&) e* Hst . 



e * HstP ( Eh ' ,1 ;h-t Hst \h'\(h'\ 

n{E v ) 



By evaluating the commutators, we get the non- 
Markovian version of our master equation for environ- 
ments with discrete spectra 



P(E b ,t) P(E b ,t) A 2 



n{E b ) 



i(E b ) t? 



+e^' T (b\B K \b')(b'\ B K i \b) S K S K (-r) 



EE /**■{ ( 34 ) 

P(E b ,t) 



n(E b ) 



-e*""' T (b\B K \b')(b'\B K ,\b)S* P{E ^ S k \-t) 

n{E v ) 

- e -^'^b\B K ,\b'){b'\B K \b)S K \-T) P{E * ,t) S K 



n{E b 



+e-^' T (b\B K >\b')(b'\B K \b) 



P{E b ,t) CK 
n{E b ) 



where uj bb i = (E b — E b >)/h are the Bohr frequencies of 
the environment. 

We assume now that the spectra of the environment is 
dense enough to be treated as quasicontinuum so that we 
can use the following equivalences 

P(E b ,t)=n(E b )tv B \b)(b\p(t) (35) 
^ P(e,t) =tv B S(e- H B )p(t) , 
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The non-Markovian version of our master equation for 
environments with a quasi-continuous spectrum is there- 
fore 



We can now rewrite Eq. (|36[) as 



P(e,t) = C s P(e,t) 



(36) 



+ e i (^>/ h n(e'){e\B K \e , ){e'\B K ,\e}S K S K '(-r)P(e,t) 
-S^>^WW){e'\B K> \e)S K P{e',t)S K '{-T) 
-e-^-^/ H n(e){e\B K ,\e')(e'\B K \e)S K '(-T)P(e',t)S K 
+e-^-^IMz')(ABAt')(^B K \e)P{e 1 t)S*\-T)S K ) . 

This equation was first obtained in Ref. Notice that 
the reduced density matrix of the subsystem is obtained 
from the quantity P(e,t) using 



Ps(t)=tv BP (t) = / deP(e,t) 



(37) 



The quantity P(e, t) can be seen as an environment en- 
ergy distributed subsystem density matrix. One should 
also realize that Eq. (|3l)|) is a closed equation for the 
P(e,i)'s but cannot be converted without approxima- 
tions into a closed equation for ps{t) using (|37[) . This is a 
first indication that Eq. (|36[) contains more information 
than what can be contained in a closed equation for ps(t). 

In the environment space, the equilibrium correlation 
function between two coupling operators B K and B K i is 

<W(<) = tr BP eq c^ HBt B K c-^ HBt B K , , (38) 

where p eq is the equilibrium density matrix of the en- 
vironment. If the environment is described by a micro- 
canonical distribution at energy e 



the correlation function reads 



(39) 



CW(M) = / &'n(e')e l(£ - £ ' )t/?l (e|i? K k'>(e'|S K 'k> -(40) 



Its Fourier transform [&(u) = ^e luJt a(t)] is 

& KK >(e,w) = hn{e + huj)(e\B K \e + huj)(e + hu\B K ,\e) (41) 

This quantity has a useful physical interpretation. 
J2kk' i s ' U P to a iac tor h 2 /(ir\ 2 ), the Fermi golden 
rule transition rate for the environment in a microcanon- 
ical distribution at energy e to absorb (resp. emit) a 
quantum of energy hu> when submitted to a periodic per- 
turbation J2 k Bk, cos u)t. Using (|4*T]) , we can easily verify 
the important property 

a KR /(e, w)n(e) = a„'„(e + fiw, — w)n(e + fiw) . (42) 



P(e,i) = C s P(e,t) 



A 2 v- '* 



fir 



do; 



{ (43) 



-e-™ T a KK ,(e,iv)S K S K (-r)P(e,*) 



-e^cw K (e, W )P(e,i)S K (-r)5 K 
+e- lWT (3 K ' K (e + fiw, -w)5 K P(e + hu, t)S K ' (-• r) 
+e <WT a KK > (e + fiw, -w)5 K '(-^)P(e + ^, t)S K \ . 

This equation is the same non-Markovian master equa- 
tion as (|36j) but written in a more compact and intu- 
itive form. It is explicit now that the effect of the en- 
vironment on the subsystem dynamics only enters the 
description via the environment microcanonical correla- 
tion function. In standard master equations, the same is 
true but with canonical instead of microcanonical corre- 
lation functions. The Markovian version of this equation 
is obtained by taking the upper bound of the integral to 
infinity J" — > f£°. This approximation is well known in 
the literature and is justified when the environment cor- 
relation function decays on time scales t c much shorter 
then the fastest time scale of the free subsystem evolution 
ts (typically given by the inverse of the largest subsystem 
Bohr frequency). 

A remark concerning the terminology is required at 
this point. In our terminology, the generator of time 
evolution is said non-Markovian if it explicitly depends 
on time whether or not it acts on ps or P(e). An alter- 
native definition (see [33j]) defines the system dynamics 
as non-Markovian if it cannot be described by a time- 
independent generator acting on p$. With this definition, 
Eq. (|4*5|) with the upper bound of the time integral taken 
to infinity would describe a non-Markovian dynamics. 



III. ROTATING WAVE APPROXIMATION 

In this section, we will use the rotating wave ap- 
proximation (RWA) which consist in time averaging 
the Markovian version of the master equation in the 
interaction form. When applying this approximation to 
Eq. (H2J), the equation takes a simple form which allows 
to prove important results. 

The RWA is most commonly used in quantum optic 
[1, [H, EH because the free subsystem dynamics gener- 
ally evolves on times scales t$ which are much faster than 
the relaxation time scales T r induced by the coupling to 
the environment. The master equation in the interacting 
representation evolves then very slowly compared to the 
Bohr frequencies of the subsystem which can therefore 
be averaged out. In other words, the RWA is justified 
if the time scale separation t c < < r r exist. In the 
mathematical physics, peoples usually refer to this aver- 
aging procedure (which is performed after a rescaling of 
time t' — \ 2 t) as the weak coupling limit [13, [H| since 
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the smaller the coupling the longer r r . Their motivation 
is to impose a Lindblad form to the master equation gen- 
erator in order to ensure the positivity of the subsystem 
density matrix [1, [ll|, |29| . The same situation occurs 
in our case. In Ref. [33|, Breuer has generalized the 
Lindblad theory to generators which act on projected 
total density matrices of the type PD|) which correlate 
the system-reservoir dynamics. Eq. (|43p is not of the 
generalized Lindblad form and could in principle lead to 
positivity breakdown similarly as the Redfield equation 
0, [H, III H|. The RWA can be used to guaranty that 
our equation preserves the positivity of the subsystem 
density matrix. By writing Eq. (|43[) in the interaction 
representation and by projecting the resulting equation 
in the subsystem eigenbasis, we get 



PA At) 



ft 2 



POO poo 

^VL dr LA 

KK' 88' 



(44) 



-e v ^ «« > e •' a KK ,{e,uj)b sg b ss ,P s , sl (e,t) 

A(u)—uj--t)TAu)- i! it i \ tjI ( i\ nft' nK 



_|_g — i('^+ u 's's') T Q i ( U 'ss+ UJ s's') t 



a KK , (e + tiu, -u)S^P*,, (e + hu, t)S?, s , } 



The RWA consist in time averaging limy—xx, ^ J^ T dt 
the right hand side of Eq. (|44|) so that 



Markovian version of our master equation in the RWA 
PUe,t) = ^EE{ ( 46 ) 

kk' ss' 

-7tS S s'& kk ' (e, -ui S s)S^Sg s P^, (e, t) 

-TrS Ss ,a K , K (e, uj s , s> )P^ (e, t)S^,S?, s , 

+7r [(1 - S ss ')S sS Ss' S ' + Sss'Ss'g] a K < K (e - fiWs's', Vs's') 

s:,piAe-fus> s >,t)s?: s , 

+TT [(1 ~ 5ss')SssSg's> + Sss'Ss's] a KK '(e + hu sS , -UJ sg ) 
SsgPss' ( £ + tu^ss, t)Sg, s , 



+iS eS > 



oj + UJ Ss 



, oj - uj s _ 
-i [(1 - 5 ss >)5 sS 5s' s < + 6 SS 'S S ' S ] Ss S 



LU + L0g>, 



+i [(1 - 5 ss i)5 sg 5 s i s > + <5 s;j '(5, 
duuP 



OK 

s'sj *J s j 



ow(e + friu, -u))P~, ( e + ^ 



This equation is a central result of this paper. It might 
look a little complicated but is in fact relatively simple. 
The first four terms are responsible for the damping of the 
subsystem and the four last term are small shifts in the 
Bohr frequencies of the subsystem. The populations (s = 
s') evolve independently from the coherences (s =/= s') 
and the coherences evolve independently for each other 
following exponentially damped oscillations 

P ss , (e, t) = (-r„/ - in ss , )P BS , (e, t) , (47) 

where the damping rates are given by 

r »' = ^E{- 27r5 -'( £ '°)^X K ' S ' (48) 



h 2 



POO poo 

^VL dr LA 

KK, SS' 



(45) 



-^e-^+^cw (e, ~)S^S:jA (e, t) 
-6 Ss ,e^-^'^a K , K (e,u;)PAAt)S*; s ,S*, s , 
+ [(1 - 5 SS ,)6 SS 6 S , S , + 5 SS ,S S , S ] e-^+*w)T 

a K > K {e + hiu, -u)Ss S PA (e + hw, t)Sg, s , 
+ [(1 - 5 ss ,)5 sS 5? s > + 5 ss ,6s, s ] e^~^ T 

a KK - (e + fej, -u)S£pL, (c + ^, t)5^ y } 



00 j„±iwt 



Using J dre 



nS(ui) ± iP — , we finally get the 



+7T 



E 



and where the modified Bohr frequencies are given by 



A 2 

J-? EE 



w + w, 
a KK '(e, w) 



s' S ss' 



.(49) 



This equation is local in the energy of the environment. 
This is not the case of the equation which rules the pop- 
ulation dynamics 

p Ae,t) = ^EE{ (so) 

kk' s 

-27ra KK ,(e,-L>s S )S«sS% s P ss (e,t) 
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This equation is a kind of Pauli equation for the total 
system which preserves the unperturbed energy of the 
total system 



(H )t = f de£(£ s +e)P Sfl (e,t) 
= dee f(e) , 



(51) 



where 



f{e)=Y,P a .{e-E a ,t) 



(52) 



represents the probability distribution inside a given en- 
ergy shell of energy e of the unperturbed total system. 
Indeed, using Eq. ([50|) , we find 



/(e) = ^P ss (e-P s ,i) = 



(53) 



This shows that the energy of the total system is con- 
served by the dynamics because the probability is pre- 
served inside each energy shell of the total system. If 
the initial condition of the total system is a product 
of a subsystem pure state of energy E s with a micro- 
canonical distribution at energy eo for the environment 
[P(e, 0) = |s)(s|(5(e — eo)] the energy distribution of the 
total system corresponds to /(e) = 8{e — E s — eo). If 
the subsystem is initially described by a density matrix 
p s (0), we get /(e) = J2 S ps,ss(0)5(s- E s - e ). The pop- 
ulation dynamics independently evolves in each energy 
shell of the total system. The equilibrium distribution of 
Eq. flU is 



P ss (e - hhj Ss , oo) 
P ss (e,oo) 



w'(e 



Using Eq. (J42J) and e = e + E s , Eq. ([54]) becomes 

P g g(e - E g , oo) _ n(e - E s ) 
P ss (e-E s ,<x) ~ n(e-E s ) ' 



(54) 



(55) 



which expresses detailed balance at equilibrium. Because 
/(e) is invariant under the dynamics, we finally get 



P ss (e -E s ,oo) 



n{e - E s ) 
E s n(e-E- S ) 



(56) 



At equilibrium, each quantum level inside a given energy 
shell of the total system has the same probability. This 
means that our equation describes how any initial dis- 
tribution inside a given energy shell of the total system 
reaches equilibrium. 



IV. EQUIVALENCE BETWEEN THE 
ENSEMBLES 



In this section, we start by showing using a simple qual- 
itative argument that if the environment is initially de- 
scribed by a canonical distribution and is assumed large 



enough for not being affected by the subsystem dynam- 
ics, our equation reduces to the Redfield master equa- 
tion. However, the most important part of this section 
is devoted to the problem of understanding in detail how 
our master equation (which rules the dynamics of a sub- 
system interacting with an environment described by a 
microcanonical distribution) effectively reduces to a Red- 
field master equation (which rules the dynamics of a sub- 
system interacting with an environment described by a 
canonical distribution) if the equivalence between the mi- 
crocanonical and the canonical ensemble is satisfied for 
the environment. 

The canonical density matrix of the environment is 
related to the microcanonical density matrix of the envi- 
ronment by 



P eq 

r can 



where 



W(e) 



de W{e) pi 



-0e 



n(e) 



(57) 



(58) 



We notice that the condition (|42|) found for the micro- 
canonical correlation functions is in fact at the origin of 
the KMS condition for the canonical correlation func- 
tions. Indeed, using (|4"2"| and 



& KK r(j3,w) = J de W(e) a K ^(e,u) , 
we obtain the KMS condition 



(59) 



(60) 



A very easy way to see the link between our equation 
and the Redfield equation is to assume that the envi- 
ronment is initially described by a canonical distribution 
and remains in it at any time during its interaction with 
the subsystem. This assumption can be expressed by the 
ansatz 



P(e,t) ^ p s (t)- 



n(e) 



Z 



(61) 



This can be qualitatively justified by assuming that the 
environment is very large compared to the subsystem. By 
integrating Eq. (|43|) over the energy of the environment 
and with the help of Eq. (ED, Eq. (|59")1 and ((BTj) . we get 
the Redfield equation [j, Hi 



p s (t) = C sP s(t) + 



A 2 ^ '* 



dr\ 



(62) 



KK 

-cw 09, r)S K S K ' (-r)ps(t) - a* KK , ((3, r)p s (t)S K ' (-r)S K 
+a* KK , ({3, T)S*p s (t)S K ' (-r) + cw (/?, r)S K ' (-r)p s (t)S K } 

witha(r) = J dioe~ iuJT a(uj) . Exactly the same procedure 
can be used for Eq. (|4"6"]) and we thus get the RWA version 
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of the Redfield equation (also called the weak-coupling- 
limit master equation) [1, EH, E2, EH • In this last case, 
one can also show that, by integrating the equilibrium 
distribution (|54[) over the energy of the environment and 
by using (fBTj) and the KMS condition ([60]) , we get 



Ps 

pZ 



(63) 



As expected, using the normalization condition £\ p§ g = 
1, we find that the subsystem equilibrium distribution of 
the RWA form of the Redfield equation is a canonical 
distribution at the temperature of the environment 



-{3E S 



Zs 



(64) 



where Zs = ^2 e~ l3E " ■ 

After this qualitative discussion, we now show that the 
precise condition for the Redfield equation to provide an 
effective description of our master equation is the equiv- 
alence between the canonical and microcanonical ensem- 
bles. Integrating Eq. (|43|) over the energy of the envi- 
ronment and using (|37p . we get 



PS 



\2 ft i-oo 

(t) = C sP s(t) + ¥ J2J Q dT J **{ ( 65 ) 



de cw(e,^)S K S K (-r)P(e,*) 



— e 



+e" 



+e" 



de ■ 



;(e,Lu)P(e,t)S K (-t)S k 



de c^ K (e,- W )S K P(e,t)S K (-r) 
de a KK ,(e,-tu)S K '(~T)P(e,t)S K } 



In order to close the equation for the subsystem density 
matrix, we have to assume that the microcanonical cor- 
relation functions can be put out of the energy integral. 
This can only be justified if the environment is such that 



a KK > (e + hw s , w a KK > (e, uj) 



(66) 



By doing this, we obtain Eq. (l62j) . but where the 
canonical correlation functions have to be replaced by 
the microcanonical correlation function a KK >(/3,u)) — > 
a KK r(e,u>). Therefore, in order to obtain the Redfield 
equation, we need to further assume that we can identify 
the microcanonical with the canonical correlation func- 
tions of the environment 



(67) 



Let us now find the conditions under which the two as- 
sumptions (|66p and (|67[) are valid. The entropy, associ- 
ated with the microcanonical distribution in an energy 
shell of the environment of width 6 e and corresponding 
to an energy e, is given by 



5(e) = k B In w — k B lnn(e)<5 £ , 



(68) 



where w is the complexion number (i.e., the number of 
available states in the energy shell). The microcanoni- 
cal temperature of the environment associated with this 
microcanonical entropy is defined as 



(3(e) 



1 



k B T(e 



1 dS(e) _ 1 dn(e) 
ks de n(e) de 



(69) 



Suppose that the environment is in a microcanonical dis- 
tribution at the energy e m . This environment can be 
effectively described by a canonical distribution at the 
temperature (fcs/3) -1 if W(e) is a sharp function with its 
maximum at the energy e m , which is therefore the most 
probable energy. In this case, we can expand hilt^e) 
around e m . We get 



InW(e) = \nW(e m ) + 
rf 3 



de 3 



de 2 
biW(e) 



lnW(e) 



(e - e m f 



2! 



3! 



,(70) 



since d\nW(e m )/de = 0. Because the energy e m corre- 
sponds to the maximum of lnW^(e), the canonical tem- 
perature is equal to the microcanonical temperature at 



- /3(Cm) • (71) 
Now using the microcanonical heat capacity defined by 

1 _ dT{e) 



(7(e) 



de 



we have that 
d 2 



de 2 



]xiW{e) 



d/3(e) 
de 



k B T 2 (e)C(e) 



(72) 



(73) 



We see now that the rational to truncate the series (|7Q|) is 
a large and positive heat capacity. This is true for most 
environments in the limit of a large number of degrees 
of freedom N, since typically C(e) ~ N. We can now 
rewrite Eq. ([70)) as 



W(e) « W{e m ) exp 



{e-e m f 



2a 2 (e m ) 



where 



c 2 (e m ) = k B T 2 (e m )C(e m ) 



(74) 



(75) 



This result confirms that e m is the maximum of W(e), 
and also shows that e m is its mean value if a(e m ) is small 
compared to the typical energies of the environment. 
This is true for large N, since e/a(e m ) ~ Vn. Under 
these conditions, an environment described by a micro- 
canonical distribution at the energy e m can be effectively 
described by a canonical distribution at the temperature 
P(e m ) so that the second assumption (|BTj) becomes justi- 
fied. However, in order to justify the first assumption 
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(l66|) . the environment effective canonical temperature 
has also to remain unchanged under energy shifts of the 
environment energy of the order of the typical subsystem 
quanta hu>s 

/?(<= ± tius) = /3(e) ± ^ hco s + ... . (76) 
The condition to have such an isothermal environment is 



frw s 


d(3{e) 






/3(e) 


de 




T(e)C(e) 



< 1 • (77) 



This condition is again satisfied when the number of de- 
grees of freedom of the environment becomes large be- 
cause hcos /(T(e)C(e)) ~ 1/N. Our conclusion is that, in 
the limit of an infinitely large environment N — > oo, cr(e) 
becomes infinitely small compared to typical environ- 
ment energy scales (so that we have the equivalence be- 
tween the microcanonical and canonical ensembles) but 
also becomes infinitely large compared to typical sub- 
system energy scales fkus (so that the environment is 
isothermal for the subsystem) . It is in this limit that the 
Rcdfield master equation becomes valid and can provide 
an effective description of our master equation. 

V. CONCLUSIONS 

In this paper, we have considered a subsystem inter- 
acting with an environment which has a sufficiently large 
number of degrees of freedom so that its spectrum can 
be supposed quasi-continuous. As a consequence, this 
environment can drive the subsystem into a relaxation 
process on time scales typically shorter than the Heisen- 
berg time of the environment [tu — Kn(e), where n(e) is 
the average density of states of the environment]. How- 
ever, the number of degrees of freedom of this environ- 
ment may still be too small for the transitions to leave 
it isothermal. Indeed, the heat capacity is finite and the 
energy exchanges between the subsystem and the envi- 
ronment can modify the energy distribution of the en- 
vironment. This is not an unrealistic assumption since 
the density of states of the environment grows exponen- 
tially with the number of degrees of freedom albeit the 
heat capacity only grows linearly. By using projection 



operators, we derived a master equation governing the 
relaxation dynamics of such a subsystem. This equation 
describes the evolution of the subsystem density matrix 
distributed over the energy of the environment. This al- 
lows us to take into account the changes of the environ- 
ment energy distribution due to its interaction with the 
subsystem. By performing the rotating wave approxima- 
tion (RWA) on this master equation, we have been able 
to show that the subsystem populations evolve indepen- 
dently from the coherences. The coherences decay in 
the form of exponentially damped oscillations, while the 
populations obey a kind of Pauli equation for the total 
system. This equation conserves the energy of the to- 
tal system and its equilibrium distribution corresponds 
to the uniform distribution of the probability over the 
energy shell of the unperturbed total system. Our equa- 
tion provides a natural representation of the dynamics of 
a subsystem interacting with an environment described 
by a microcanonical distribution. The Rcdfield master 
equation is the usual way to represent the dynamics of a 
subsystem interacting with an environment described by 
a canonical distribution and is a closed equation for the 
density matrix of the subsystem. We have shown that, if 
the equivalence between the microcanonical and canoni- 
cal ensembles is satisfied for the environment (containing 
infinitely many degrees of freedom) , our equation reduces 
to the Redfield master equation. If this equivalence is not 
satisfied, our equation becomes of crucial importance to 
correctly describe the subsystem relaxation. In the same 
sense that the microcanonical ensemble is more funda- 
mental than the canonical ensemble, our master equa- 
tion is more fundamental then the Redfield equations. 
We believe that this work improves the understanding 
of kinetic processes in nanosystems where the thermody- 
namical limit cannot always be taken. 
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